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CN ■ The one parameter scaling theory is a powerful tool to investigate Anderson localization effects 

in disordered systems. In this paper we show this theory can be adapted to the context of quantum 
, chaos provided that the classical phase space is homogeneous, not mixed. The localization problem 

' in this case is defined in momentum, not in real space. We then employ the one parameter scaling 

theory to: a) propose a precise characterization of the type of classical dynamics related to the 
Wigner-Dyson and Poisson statistics which also predicts in what situations Anderson localization 
corrections invalidate the relation between classical chaos and random matrix theory encoded in 
the Bohigas-Schmit-Giannoni conjecture, b) to identify the universality class associated with the 
, metal-insulator transition in quantum chaos. In low dimensions it is characterized by classical 

^ ■ superdiffusion, in higher dimensions it has in general a quantum origin as in the case of disordered 

systems. We illustrate these two cases by studying Id kicked rotors with non analytical potentials 
and a 3d kicked rotor with a smooth potential. 
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PACS numbers: 72.15.Rn, 71.30.+h, 05.45.Df, 05.40.-a 

d ■ 

[ The Bohigas-Giannoni-Schmit (BGS) conjecture [l| states that quantum spectra of classically chaotic systems are 

^ universally correlated according to the Wigner-Dyson random matrix ensembles (WD). Therefore for scales of time 

. comparable to the Heisenberg time the quantum system evolves to a state where all dynamical features are washed 

^ ■ out, and all the available phase space is occupied. Eigenstates are fully delocalized and the spectral correlations are 

^ ' dictated by symmetry only. For the sake of completeness we mention the relation between level repulsion and classical 

lO chaos was noted previously in [2, 3]. However we stick to the term "BGS conjecture" through the paper since O. 
OO ■ . 

■ Bohigas and coworkers were the first to claim that a) this relation was universal, b) spectral correlations are described 
by random matrix theory. From our point of view this is a qualitative difference from the findings of Ref. 

SB. 

^ ; n, 

I Similarly, according to the Berry- Tabor-Gutzwiller (BTG) conjecture [J], the spectral correlations of classically 

, integrable systems are well described by Poisson statistics. According to the Nocther theorem there is a conserved 

' quantity (generalized momentum) for each classical symmetry. Integrability assures that the number of these 

' conserved momenta is equal to the dimensionality of the space. In quantum mechanics each conserved quantity 



leads to a good quantum number. A given quantum state is thus characterized by a unique set of good momentum 
quantum numbers. Since different states are characterized by a different set of quantum numbers, the spectrum 
is not correlated. In other words, in integrable systems the Hamiltonian is diagonal in a basis of momentum eigenstates. 



The influence of these two conjectures in the establishment and development of quantum chaos is hard to 
overestimate. The reason is simple: Comparatively it is straightforward to determine whether the classical dynamics 
is chaotic or integrable. With this input one can predict the quantum transport properties of the system even without 
a full consideration of the specific details of the Hamiltonian. 
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However, exceptions are well documented in the literature, among others we cite arithmetic billiards [5|, Harper 
model and the ubiquitous kicked rotor Q], namely, a free particle in a circle which experiences periodic kicks of 
amplitude given by a periodic smooth potential, 

2 

n = ^-Kcos{q)Y,S{t-Tn). (1) 

n 

For short time scales, both quantum and classical motion is diffusive in momentum space. However, quantum 
diffusion is eventually arrested due to destructive interference effects that localize eigenstates in momentum space. 
In this limit spectral correlations are described by Poisson but not by WD statistics. This counter-intuitive feature, 
usually referred to as dynamical localization [2], was fully understood fS] after mapping the kicked rotor problem 
onto a short range one dimensional disordered system where localization is well established. 

Deviations from the BGS conjecture are also expected [Q] for eigenvalue separations 6E ^ Ti/tE due to weak 
localization effects. The typical scale = A~-'^| log?i| is the Ehrenfest time with A the classical Lyapunov exponent. 
These corrections do not invalidate the BGS conjecture since 5E ^ Ti/Ie is much larger than the mean lever spacing 
A, the natural scale in which WD statistics holds. 

Another potential source of deviations from WD statistics is related to energy scales larger than fi/tfi where tfi is 
the typical time for the classical dynamic to become fully chaotic. In billiards it is just the flight time of a particle 

10, ll| in terms of the spectral determinant of 



across the system. These corrections can be expressed analytically 
the Perron-Frobenious operator which controls the relaxation of the classical density to the ergodic-universal limit 
described by random matrix theory. As for weak localization corrections, the associated energy scale is much larger 
than the mean level spacing and consequently does not really affect the validity of the BGS conjecture. 

From the previous discussion it is clear that the main deviation from the semiclassical picture contained in the 
BGS conjecture is due to strong quantum localization effects, or as it is usually referred in the literature, Anderson 
localization. After this introduction some questions arise: Why in some quantum chaotic systems localization is 
only weak and in others is strong enough to make the system and insulator invalidating thus the BGS conjecture? 
Are chaotic systems the only ones in which the BGS conjecture applies? If not so, is it possible to give a precise 
relation between classical motion and quantum features? And is it possible to define any other universality class in 
the context of non random Hamiltonians? 



In this paper we aim to address these questions by using ideas and techniques originally developed in the theory 



of disordered systems [12[ such as the one parameter scaling theory [l3j. 



We start with a brief review that highlights some of the main ideas and results of this field. One the main problems 
in the theory of disordered systems is the study of the role of quantum effects such as tunneling and interference on 
the transport properties of a particle in a random potential. In general, transport properties are strongly affected by 
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both the dimensionahty of the space ,d, and the strength of disorder. 



In two and lower dimensions destructive interference caused by backscattering 



produces exponential localization 



of the eigenstates in real space no matter the amount of disordered considered |12l |. As a consequence, quantum 
transport is suppressed, the spectrum is uncorrelated (Poisson) and the system becomes an insulator. In three and 
higher dimensions there exists a metal insulator transition for a critical amount of disorder. Thus for disorder below 
the critical one the wavefunctions are extended through the sample, namely, they are effectively represented by a 
superposition of plane waves with random phases. The Hamiltonian is accurately approximated by a random matrix 
with the appropriate symmetry, and the spectral correlations are given by WD statistics In the opposite limit, 
wave functions are exponentially localized as in d < 2, and Poisson statistics applies. 

As in the quantum chaos case these universal features are restricted to long time scales (related to energy scales 
of the order of the mean level spacing) such that an initially localized wave-packet has already explored the whole 
phase space available. In other words, universality is related to certain ergodic limit of the quantum dynamics. For 
shorter time scales the system has not yet relaxed to the ergodic limit and deviations from universality are expected. 



3. 



These features can be understood within the framework of the one parameter scaling theory [13] . A key concept 
in this theory is the dimensionless conductance g introduced by Thouless [isj l. It is defined either as the sensitivity 
of a given quantum spectrum to a change of boundary conditions in units of the mean level spacing A cx or as 
g = Ec/ A where, Ec, the Thouless energy, is an energy scale related to the diffusion time to cross the sample. In 
the semiclassical limit Ec = fiDda/ {Dda is the classical diffusion constant) and therefore g oc U^^^ . On the other 
hand if the the particle is exponentially localized, g oc er^l^ where ^ is the localization length and L is the system 
linear size. 

The change of the dimensionless conductance with the system size is an indicator of localization. Qualitatively, in 
the L — > cx) limit, g — > 00 in a metal and g — > in an insulator. In order to proceed it is useful to define 

«»)^^ (^) 

which describes the running of g with the system size. 

From the above definitions, fi{g) = (i — 2>0ina metal (without quantum corrections) and fi{g) = \og{g) < in 
an insulator. 

The one parameter scaling is based on the following two simple assumptions: 

• The l3{g) function is continuous and monotonous; 

• The change in the conductance with the system size only depends on the conductance itself which is a function 
of the system size. 



We are now ready to analyze how localization corrections depend of the spatial dimensionality d: 
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1. d =1. For g — > oo, f3{g) = — 1 < 0. In the opposite limit g — > 0, (3{g) = log{g) < 0. Since f3{g) < for any g, 
the dimensionless conductance always decreases with the system size and the system will be an insulator in the 
L — !■ oo no matter the amount of disorder. 

2. d =3. For g — > oo, I3{g) = 1 > 0. In the opposite limit however g 0, /3(g) ~ log(.g) < 0. Since j3{g) 
is continuous and monotonous there must he g — gc such that P{gc) — 0. At g ~ gc the system undergoes 
a metal-insulator transition. The one parameter scaling theory thus predicts a new window of universality 
characterized by a scale invariant dimensionless conductance g ~ gc- For a fixed disorder, if g{L) > {<)gc then 
the system will flow to the metallic (insulator) limit as i — > oo. A similar picture holds for c? > 3. The spectral 
correlations at the Anderson transition, usually referred to as critical statistics [igI . 17 |, are scale invariant and 
intermediate between the prediction for a metal and for an insulator [itI . fisl ] . By scale invariant we mean any 
spectral correlator utilized to describe the spectral properties of the disordered Hamiltonian does not depend on 



the system size. Eigenfunctions at the Anderson transition are multifractals Q, (for a review see [2ll. l22j|). 
namely, their moments present an anomalous scaling, Vq ~ f d'^r\tp{r)\'^'' oc L^^i'-i^^'' with respect to the 
sample size L, where Dq is a set of exponents describing the Anderson transition. 

3. d —2. For g ^ oo, (3{g) = but for g — > 0, [3{g) = log(g) < 0. In this case, since f3{g) = 0, the semiclassical 
analysis is not enough to determine the localization properties in the thermodynamic limit. We have to determine 
whether (3{g) is positive or negative (or it remains zero) once quantum interference effects (weak localization 
corrections) are included. Technically one uses standard perturbation theory to compute one loop quantum 
corrections to the classical diffusion constant included in the definition of the dimensionless conductance. The 
result of this calculation is that I3{g) ~ —a/g {a > 0) for g ^ 1. Therefore P{g) < for all g and, as in the 
Id case, all eigenstates are localized in the L oo for any amount of disorder. We note that the negative sign 
of the weak localization correction is not universal, it depends on the details of the Hamiltonian. We shall see 
that if the classical diffusion is anomalous, weak localization corrections renormalize the value of the classical 
diffusion constant, but the semiclassical critical point f3{g) = is not modified. 

A few comments are in order: 

1. The Anderson transition in a disordered conductor is a consequence of a highly non trivial interplay between 
quantum destructive interference effects and quantum tunneling. In low dimensions, d ^ 2, weak quantum 
destructive interference effects induce the Anderson transition. Analytical results are available based on per- 
turbation theory around the metallic state In high dimensions, d 3> 2, quantum tunneling is dominant 
and the locator expansion or the the self-consistent formalism can be utilized to describe the transition. 
We note that in these papers corrections due to interference of different paths are neglected. 

2. The above analysis is strictly valid only for time reversal invariant systems with no interactions. A strong 
magnetic field or a spin-orbit interaction may induce a metal insulator transition in two dimensions. 

To summarize: According to the BGS conjecture, the level statistics of quantum systems whose classical counterpart 
is chaotic are described by WD statistics and eigenfunction are delocalized in momentum space. According to the 
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FIG. 1: (Leh)P{p,t) versus p (log scale). In all cases t = 2000 except for a = -0.5 where t = 40. For V{q) = with a = 0.9 
(solid) and V(q) = {\q\ + 0.005)"'^ (dash) the diffusion is normal (squares and circles are for best Gaussian fitting respectively). 
However for a = 0.4 (dot) and a = —0.5 (dash-dot), P{p,t) oc tp~'^°'\ Inset shows 7 as a function of a. The star corresponds 
to V{q) = log(|g|). (Right) Moments of the classical distribution for different a's. As a increases classical diffusion is slower. 



BTG conjecture, the level statistics of quantum systems whose classical counterpart is integrable are described by 
Poisson statistics and the eigenstates are localized due to symmetry constraints. Strong localization effects typically 
invalidate the relation between classical chaos and WD statistics. 

In the context of disordered systems the one parameter scaling theory provides a valuable framework to understand 
localization effects. It predicts three universality class: For g ^ 00 the system is a metal, eigenfunctions are delocalized 
in real space and the spectrum is correlated according to WD statistics. For g ^ the system is an insulator, 
eigenfunctions are exponentially localized and the spectrum is correlated according to Poisson statistics. For g — gc 
a metal insulator-transition takes places. The eigenstates are multifractal and the spectral correlations are universal 
(critical statistics) but different from WD and Poisson statistics. 

The organization of the paper is as follows: In section 2 we discuss under what conditions the one parameter scaling 
theory can be applied to quantum chaos. Then we identify a variety of problems in quantum chaos which the one 
parameter scaling theory may help solve. In section 3 we tackle one of them: the universality class related to the 
Anderson transition in quantum chaos. We test numerically the predictions of the one parameter scaling theory about 
the Anderson transition in quantum chaos in a 3c? kicked rotor with a smooth potentials and in a Id kicked rotor with 
singularities. Section 4 is devoted to conclusions. 



I. ONE PARAMETER SCALING THEORY IN QUANTUM CHAOS 



In this section we explain how to adapt the one parameter scaling theory to the context of quantum chaos. 
The absence of a true ensemble average in quantum chaos is in principle a major obstacle for a straightforward 
application of techniques and ideas originally developed in the field of disordered systems. Fortunately the one 
parameter scaling theory is an exception. The basic assumptions (continuity and monotonicity of the f3{g) function 
and its sole dependence on g) do not lie on any ensemble average. Moreover the mean level spacing in the 
definition of the dimcnsionlcss conductance have a well defined meaning without having to perform an ensemble 
average. The Thouless energy, which is related to the classical diffusion time through the sample, is typically 
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FIG. 2: Comparison of the classical and quantum time dependence of the moments {\p\'") oc t'^''^^ for the case of a logarithmic 
singularity (a = 0) related to the metal insulator transition. The classical time dependence is not modified by quantum 
corrections (right). However the value of the diffusion constant is smaller in the quantum case (left). This is an indication 
that quantum interference effects slow down the classical motion but cannot induce full dynamical localization. The same 
conclusions apply for a < 0. 

understood as an average over many realizations. However this is not necessary. We just need to find a way 
to determine a typical and meaningful travel time through the sample. For instance, in kicked system it is the 
time it takes to explore a momentum basis of a given size. The main problem in quantum chaos is that in mixed 
systems it will depend strongly on the initial conditions. If the particle is initially located in a integrable region 
the typical time will be larger than in a chaotic region since KAM barriers slowdown transport in momentum 
space. Thus in order to define properly the dimensional conductance in quantum chaos we have to impose the re- 
striction that the classical phase space must be homogeneous. In this way a travel time can be defined unambiguously. 

Another issue that requires clarification is the space in which the localization problem is defined. In disordered 
systems is always the real space. Typically one fixes the strength of the disordered potential (for instance the density 
of scatterers) and then study how the dimensional conductance vary with the system size. By contrast, in quantum 
chaos the localization problem is typically defined in momentum space. Therefore the Thouless time is defined as 
the time that a wavepacket needs to explore a finite momentum basis of size N (equivalently to the system size L in 
disordered systems). As in disordered systems this can be estimated by computing the second moment {p^) of the 
distribution. 

The reason why localization is naturally defined in momentum space is related with the Noether theorem in 
classical mechanics which states that for each classical symmetry there is a conserved quantity. If a system is 
classically integrable the number of conserved quantities (canonical momenta) is equal to the dimensionality of the 
system. Semiclassically it is evident that each of these canonical momenta becomes a good quantum number which 
labels the state. An integrable system is localized in momentum space since there exists a basis of good (momentum) 
quantum numbers in which the Hamiltonian is diagonal and consequently the spectral correlations are Poisson-like 
as for an insulator. If the number of conserved momenta is less than the dimensionality of the space, the matrix 
representation of the Hamiltonian in any basis of momenta become non-diagonal and the spectrum is correlated to 
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some degree. According to the BGS conjecture, in the case of no symmetry the spectral correlations are sufficiently 
strong that WD statistics applies. We note that this semiclassical picture does not take into account of strong 
localization effects that tend again to drive the system to the insulator limit. 

To summarize: in quantum chaos the dimensionless conductance is a meaningful quantity provided that clas- 
sical phase space is homogeneous. The localization problem in quantum chaos is naturally defined in momentum space. 

The classical dynamics in quantum chaos is by no means restricted to standard diffusion (p^) = Dt. Different types 
of anomalous diffusion, 

- (3) 

with /? 7^ 1 may be relevant. The mean level spacing is in many cases given by A 1/N'^, but there are important 
exceptions: i) for periodic potentials, the Bloch theorem applies, the spectrum is continuous, and A = 0; ii) for 
systems whose eigenstates are exponentially localized, A ^ even in the N oo limit; iii) in systems with a singular 
continuous spectrum the scaling with the system size may be anomalous 

A cx 7V-d/de (4) 

(de 7^ 1). A precise definition of de may depend on the system in question. In the Harper model c?e ~ 1/2 stands 
for the Hausdorff dimension of the spectrum. 

We are now ready to define the dimensionless conductance in quantum chaos. In cases i) and ii) above, g —^ oo 
(metal) and g — (insulator) respectively. In case iii) (including de — 1), we define, with the help of Eq. ([3]) and ([4]), 
the dimensionless conductance in quantum chaos. 



In disordered systems it is intuitively clear that if the strength of disorder is sufficiently strong the insulator limit 
will be eventually reached. In quantum chaos it is not that evident as it may occur that the coupling constant of the 
model which plays the role of disorder never reaches this limit. This is a fact to take into account before making any 
statement about the existence of a metal-insulator transition in quantum chaos. 

In order to study the change of the dimensionless conductance with the system size due to quantum effects, it is 
helpful to define the function l3{g): 

m - (6) 

ologL 



If quantum effects are not important. 



Pig) = Idas = X ~ ^- 

de P 
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Based on the above expressions and under the assumptions of the one parameter scaUng theory we propose the 
following alternative definition of universality class in quantum chaos: 

For systems such that jcias > 0, eigenfunctions are delocalized as in a metal and the spectral correlations are 
universally described by WD statistics. 

For systems such that jcias < 0, eigenfunctions are localized as in a insulator and the spectral correlations are 
universally described by Poisson statistics. 

For systems such that jcdas = 0, eigenfunctions are multifractal as at the metal insulator transition and the spectral 
correlations are universally described by critical statistics jl6| . 

Assuming for the moment that the spectrum is not fractal, de = 1, the one parameter scaling theory predicts that 
if the classical motion is such that (3 < 2/d the quantum properties are those of an insulator no matter whether the 
classical dynamics is chaotic or not. For instance, for d ~ 3, any classical motion leading to /3 < 2/3, even if chaotic, 
will drive the quantum system to the insulator phase. 

Another example is the Id kicked rotor with a smooth potential [3], classically (p^) cx t ^ "fdas = — 1 < 0. The one 
parameter scaling predicts that the spectral correlations are described by Poisson statistics as in an insulator even 
though the classical dynamics is chaotic. It is remarkable that by using scaling arguments dynamical localization can 
be predicted without having to map the problem onto a Id Anderson model. 

Equation ^ also suggests that a metal insulator transition jcias = is possible in less than two dimensions. In 
Id, it is related to a superdiffusive classical motion such that /? = 2/de > 2. For instance in non conservative Id 
systems such as the kicked rotor (c?e = 1), the spectral correlations will be metal- like for (3 > 2, similar to those of a 
metal-insulator transition for /? — 2 and, insulator-like for (3 < 2 . Therefore dynamical localization can be overcome, 
even in one dimension, if the classical diffusion is fast enough f3 > 2. 

We note that these conclusions are still restricted to the semiclassical approximation. Quantum effects such as 
destructive interference may slow down or even stop the quantum motion. However the one parameter scaling 
theory, as is the case of disordered systems, is capable to handle this situation as well. The way to proceed is to 
compute the Thouless energy not with the classical but with the quantum travel time to cross the sample. Assuming 
{p'^)quan oc t^i , WD (Poisson) Statistics applies if — d/de — 2/ [3q > (<)0. As an example we discuss the Harper 
model 6,] . It is well known that this model undergoes a metal insulator transition for a specific value of the coupling 
constant (A = 2). We note that classically the system is integrable. The semiclassical dimensionless conductance 
cannot predict the transition. However the quantum motion (x^) ^ t^"^" at the metal-insulator transition is 
diffusive where d// ~ 0.5 is the Hausdorff spectral dimension. According to our definition Eq. ([5]), the "quantum" 
dimensionless conductance g = gc is a constant independent of the system size. This is a signature of a metal 
insulator transition. Therefore our simple method predicts correctly the metal insulator transition in this model as 
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FIG. 3: The level spacing distribution P{s) in units of the mean level spacing s — {Ei — Ei-i)/A. A transition is observed 
from Poisson to WD statistics (keeping e constant) as a goes from negative to positive. In both cases = 3100. 




FIG. 4: The spectral rigidity A3(L) versus L for V{q) — elog(g). The right (left) panel corresponds to the case of (broken) 
time reversal invariance. The linear behavior for L ^ 1 is a signature of the Anderson transition. The slashed lines in the case 
of broken time reversal invariance are the analytical prediction e ~ l/ft for this type of classical anomalous diffusion. 



well. 
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A. Stability of the semiclassical predictions 

From a more technical point of view the problem of the stability of the semiclassical predictions under quantum 
effects boils down to the way in which the classical diffusion constant included in the definition of the dimcnsionlcss 
conductance gets renormalized by quantum corrections. For normal diffusion it is well known that quantum 
corrections (interference) may modify substantially the semiclassical results. For instance in a 2d disordered system 
[isi [2^ , a simple calculation of the one loop pcrturbativc quantum correction shows that the diffusion constant is 
smaller than in the classical case, I3{g) becomes negative and the semiclassical metal-insulator transition is destroyed 
by quantum corrections. 

However in the case of classical anomalous diffusion we shall show that in certain circumstances semiclassical 
predictions are stable under quantum corrections. For the sake of simplicity we focus on the case d = \ and d^ = 1. 
Extension of the results to other dimensions is straightforward. Our starting point is a classical diffusion process 
described by the following fractional Fokker-Planck equation, 

d (92//3 \ 



dt ^'"dlpiyf^ 

(see [sgI for the definition of fractional derivative) where Dda is the classical diffusion coefficient and /3 is a real 
number. The moments of the distribution P{p,t) (if well defined) are given by (IpI*^) oc t'^*''/^, with k a real numbers. 
In analogy with the case of normal diffusion the propagator is given by KQ{q,LU = 0) = — 2/0 with v the spectral 
density. 



The one-loop correction to the classical propagator due to interference effects takes the following form [20|, |22| , 

|2//3 _ 



The integral over k is carried out, with the help of the following expansion 

2qk 1 g2 2 / I \ /qk^ ^ 



22| 



q'» k 



\k\2/l3 ] \q\V0 

\k\V0 

where L, the system size, acts as an infrared cutoff. 

The quantum diffusion constant Dqua to one loop is easily obtained from Eq. ([9| , ([S]) . 

(10) 
(11) 

(12) 

where C is a different constant for each case. 

The importance of the quantum effects depends strongly on the value of 2/f3. In the region 1 < 2//3 < 2 quantum 
corrections diminish the value of the classical diffusion constant. These corrections, due to its dependence on the 



^qua 


= Dcla 


- CL^^"-' 


1< 2//3 < 2 


^qua 


— Dcla 


-Clog(<7i) 


2//3=l 


Dqua 


— Dcla 




< 2//3 < 1 
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system size, grow as higher order in perturbation theory are taken into account and eventually induce a transition to 



localization 



22| 



In the region < 2//3 < 1 the correction is of higher order in q than the classical term. This feature is not modified 

! I 

by higher orders in perturbation theory 22] and consequently the classical diffusion is not normalized. 

The case 2//3 = 1, related to the metal-insulator transition, deserves a special consideration. The logarithmic 
behavior resembles superficially that of 2d disordered system, 

Dgua= D,ia-CHL/l) , (13) 

However there are essential differences. Comparing the two formulas, we see that in the case of anomalous diffusion 
the mean free path I in the 2c? case is replaced by the inverse momentum q^^. Therefore, the correction to the bare 
coupling constant, unlike the 2D disordered conductor, is small for low momenta q ^ 1/L. Qualitatively this implies 
the absence of eigenstate localization. For a rigorous proof based on the evaluation of P{g) including higher order 
terms we refer to [22| and reference therein. 

Therefore a Hamiltonian with classical dynamics such that 2/(3 = 1 should have quantum properties similar to 
that of a disordered conductor at the metal- insulator transition. If the classical diffusion is slower (2/(3 < 1) quantum 
corrections will eventually induce a transition to localization. A final remark is in order. We recall that by dimension 
we do not mean the spatial dimension but the effective dimension in momentum space. A conservative one dimensional 
(in space) system is always integrable so no transition to delocalization can be observed. However in momentum 
space the effective dimension is zero, just a point defined by energy conservation. But if energy is not conserved, as is 
the case in kicked rotors, the effective dimension in momentum space is the same that the spatial one, namely, the unity. 



To summarize: We have identified three universality class in quantum chaos: metal, insulator and metal-insulator 
transition. According to the one parameter scaling theory these are the only windows of universality. If the classical 
diffusion is anomalous wc have shown in what circumstances semiclassical predictions are not substantially modified 
by quantum effects. Two ways to reach the metal insulator transition in quantum chaos are identified: By classical 
anomalous diffusion such that jcias = or by quantum effects 7^ in a classical Hamiltonian with normal diffusion 
and "fcias > 0. 

II. SCALING AND THE ANDERSON TRANSITION IN QUANTUM CHAOS 

We now test the predictions of the previous section in two specific models: A 3d kicked rotor and a Id kicked rotor 
with a classical singularity. In the former the classical diffusion is normal. In the latter it is anomalous with an 
exponent which depends on the type of singularity. We shall show that by tuning the kicking strength (in 3c?) or the 
type of singularity (Id) to a specific value we can induce a metal insulator transition in quantum chaos. However the 
physical origin is slightly different in both situations. In 3c? it is caused by quantum destructive interference effects 
that slow down the classical chaotic motion while in Ic? has its origin in the underlying classical anomalous diffusion. 
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FIG. 5: Classical and quantum (p^) for three different kicking values, above, below and at the transition fc = fcc ~ 3.3 in the 3d 
kicked rotor with a smooth potential. Classical diffusion is normal however quantum destructive interference slow down and 
eventually arrest the diffusion for k < kc- In agreement with the one parameter scaling theory, at k = kc, (p^) (x t'^''^ . 

Finally we mention other systems that, according to the classification scheme introduced in the previous section also 
belong to the metal insulator universality class. 



A. Anderson transition induced by classical anomalous diffusion 



Previously we determine by using the one parameter scaling theory that classical superdifFusion may induce a 
metal-insulator transition even in Id. In this section we provide numerical support to this theoretical prediction. In 
order to proceed we need to choose a system in which classical diffusion is anomalous. A suitable candidate is the 
kicked rotor with a non analytical potential introduced in [2^, 

p2 



n = ^ + V{q)J2S{t~nT) (14) 



with Vjq) = e 



see 



28 



29 



g|" and V{q) — elogd^j); q G [— tt, tt), a G [—1, 1] and e a real number (for other type of singularities 
30|). The dynamical evolution over T is dictated by the map: Pn+i ~ Pn ^ ^ a^'^"^ ' 'Zn+i ~ qn + 2^Pn+i 
(mod 27r). 

We have evaluated the classical density of probability P(p, t) and its moments {p^) by evolving 10^ initial conditions 
uniformly distributed in {q,p) G [— tt, tt) x [— 7r,7r). For a > 0.5 (see Fig. 1), the diffusion is normal and P{p,t) ^ 
g-cp I't l^/l (c a constant). However, for —0.5 < a < 0.5 and t > p, P{p,t) ^ tp'^^"-^ has power-law tails (see inset in 
Fig. 1 for a relation between 7 and a). Such tails are considered a signature of anomalous diffusion. As a consequence, 
the classical moments of the distribution scale as {p^) ^ ^(i-a)fe (ggg pjg^ and the dimensionless conductance is 



given by 



Pig) = ■ (15) 

1 — a 



According to the findings of the previous section, for a < quantum corrections do not modify qualitatively the 
semiclassical prediction for f3{g). Therefore according to the one parameter scaling theory we expect WD statistics 
for a < 0, Poisson statistics for a > (quantum corrections will only make the transition to localization faster) and 
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FIG. 6: The spectral rigidity A3{L) versus L and level spacing P{s) for three different kicking values, above, below and at 
the transition k = kc ~ 3.3 in the 3d kicked rotor with a smooth potential. A transition from WD to Poisson statistics is 
clearly observed as the kicking strength crosses the critical value k ft: kc from above. For k ^ kc the level statistics have all the 
signatures of an metal- insulator transition such as linear A3(L), level repulsion for s — > as in a metal but exponential decay 
of -P(s) for s > 1 as in an insulator. Although it is not shown the spectrum is scale invariant at fc « fcc. 



an Anderson transition for a = that corresponds naturally with the case of a log singularity. 



We now compare these theoretical predictions with numerical results. In order to proceed we evaluate the quan- 
tum evolution operator hi over a period T. After a period T, an initial state ^/^o evolves to tp{T) = Uipo = 
e ih e h e "iAo where p and q stand for the usual momentum and position operator. The spectrum was 
obtained by solving the eigenvalue problem U'i/n ~ e^"'"/''^„ where is an eigcnstate of U with eigenvalue k„. We 
express the evolution operator in a matrix form {m\U\n) in the basis of the momentum eigenstates = '^Tf^}- "^^^ 
resulting evolution matrix (for N odd) then reads 

{m\U\n) = le-'2^^'^"'/^ ^ e'-^C'™'") (16) 

I 

where (t){l,m,n) = 2tt{1 + eo)im - n)/N - iV{2T:{l + eo)/N), I = -{N - 1)/2,...{N - l)/2 and < 6*0 < 1; 6lo is a 
parameter depending on the boundary conditions (^o = for periodic boundary conditions). The eigenvalues and 
eigenvectors of U can now be computed by using standard diagonalization techniques. 

We first analyze the moments of the quantum density of of probability, (p^)- In Fig. 2 (right) it is observed 
that in agreement with the perturbative analysis above the classical time dependence is not modified by quantum 
corrections. However the diffusion constant is smaller in the quantum case due to quantum destructive interference 
effects (Fig. 2 left). Similar results are obtained for a < 0. For a > quantum corrections can only make 
the transition to localization faster. In conclusion, for a < the scaling of the dimensionless conductance does 
not change with respect to the semiclassical limit g cx L^" so a genuine metal-insulator transition is expected at a = 0. 



The analysis of the level statistics confirms this theoretical prediction. The level spacing distribution P{s) (see 
Fig. 3) tends to Poisson (WD) for a > (a < 0). As a approaches zero it is much more difficult to observe the 
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transition to either WD or Poisson statistics. The reason for that is simple. The dimensionless conductance goes as 
g{L) cx L^" but the universal WD (Poisson) limit are only expected for 5 — > oo (g ^ 0). Therefore, as we approach 
the transition a = 0, the universal limits are reached very slowly. 



We now study the case of a potential with log singularity which according to our prediction it is related to the 
Anderson transition. In order to proceed we also evaluate long range correlators as the spectral rigidity A3(_L) = 
-p- /(f (i^ - 2L'^x + x^)^/2'^{x)dx where ^^(L) = (L^) - (L)^ = L + 2 J^{L - s)R2{s)ds is the number variance and 
i?2(s) is the two level spectral correlations function. The i?2(s) associated to this type of motion has been evaluated 
in Ref.jsiJ for systems with broken time reversal invariance, 

„ , , -7 , . TT^h^ sin^fTTs) 

i?2 s = s = ^-^ — '- (17 

^ ' ^ ' 4 sinh2(7r2/is/2) ^ ^ 

where /i <C 1 is related to e by /i 1/e ^ 1. In Fig. 4 we show A3(i) for different e's and A'^ = 3099. The time reversal 
invariance in Ea. p6|) is broken by setting T = 2n(3 with /3 an irrational number. As observed in Fig. 4 (left), A3(L) 
is asymptotically linear, this feature is typical of a disordered conductor at the metal-insulator transition. Moreover, 
the agreement with the analytical prediction based on Eq. p7l) is excellent. We remark that the value of h best fitting 
the numerical result is within 5% of the analytical estimate h = 1/e. We have repeated the calculation keeping the 
time reversal invariance {T = 27tM/N) of Eq. (fT7| (Fig. 4, right). We do not have in this case an analytical result to 
compare with, but A3(L) is also asymptotically linear as in the previous case. 

These numerical results clearly support the analytical results of the previous section. Three and only three univer- 
sality class are observed: metal, insulator and metal-insulator transition. Finally we note that the novel semiclassical 
characterization of the metal-insulator transition proposed in this paper may be of interest for experimental studies 
about Anderson localization, as it enlarges the number of physical systems in which this phenomenon can be studied. 



B. Anderson transition induced by quantum fluctuations 



We now study the other route to the Anderson transition, namely, through quantum destructive interference effects 
that slow down an otherwise classically chaotic motion. 

Our starting point is again the kicked rotor introduced in the previous section. However this time it is defined in 3d 
and with a smooth potential V{9i, 62, Os) — kcos{9i) cos(6'2) cos(6'3), the classical diffusion in momentum space is thus 
expected to be normal provided that the classical space is fully chaotic. We note if the spectrum is not fractal {d^ = 1) 
and the diffusion is normal, the dimensionality must be larger than two in order to observe a metal insulator transition. 



In Fig. 5 (left) it is shown that as was expected classical motion is diffusive (p^) (x t in the range of k of interest. 
Quantum dynamics depends strongly on fc. In analogy with a 3d disordered system we expect that strong Anderson 
localization effects stop the classical diffusion for sufficiently small fc. In the opposite limit quantum effect should 
not modified substantially the classical motion. We carry out [32l | a careful finite size scaling analysis in order to 



determine how localization corrections depend on fc. A metal- insulator transition of purely quantum origin is located 
at k — kc ^ 3.3. Quantum diffusion becomes anomalous with f3 — 2/3 (see Fig. 5 (right)) in agreement with 
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the prediction of the one parameter scahng theory. In consequence the level statistics should be described by WD 
(Poisson) statistics in the limits k ^ (^)fcc. 

As is shown in Fig. 6, the numerical results fully agree with these theoretical predictions. For k = kc the spectral 
correlations have all the features of an 3d Anderson transition including the numerical value of parameters such as 
the slope of the number variance ^ 0.2. 

The kicked rotors studied in this talk are not the only ones that belong to the universality class of the metal-insulator 
transition. Other models that falls in it, at least for some value of the coupling constant, include, Harper model Q], 
quantum exchange interval maps 33|, Coulomb billiards 3J], the Kepler problem 35| and generalized kicked rotors 
28, 3- 



III. CONCLUSIONS 



The main findings of this paper can be summarized as follows: For the first time we have adapted the one parameter 
scaling theory to the context of quantum chaos. Then we have utilized it to determine the number of universality 
class in quantum chaos and propose a more accurate definition of them. The one parameter scaling theory has 
permitted us to give for the first time a precise definition of the universality class related to the metal insulator 
transition. We have identified two routes (semiclassical and purely quantum) to reach a metal-insulator transition in 
quantum chaos. In order to support our theoretical conclusions we have investigated numerically a Id kicked rotor 
with singularities where the Anderson transition has a semiclassical origin related to classical superdiffusion and the 
3d kicked rotor with a smooth potential where the transition to localization is induced by destructive interference 
effects. Our findings open the possibility of studying the metal insulator transition experimentally in a much broader 
type of systems. 



AMG acknowledges financial support from a Marie Curie Outgoing Action, contract MOIF-CT-2005-007300. JW 
is grateful to professor C.-H. Lai for his encouragement and support, and acknowledges support from Defence Science 
and Technology Agency (DSTA) of Singapore under agreement of POD0613356. 



[1] O. Bohigas, M.J. Gianonni and C. Schmit, Phys. Rev. Lett. 52, 1 (1984). 

[2] S. W. McDonald and A. N. Kaufman, Phys. Rev. Lett. 42, 1189 (1979). 

[3] G. Casati, F. Valz-Gris, and I. Guarneri, Lett. Nuovo Cimento 28, 279 (1980). 

[4] Chaos in Classical and Quantum Mechanics, M. C. Gutzwiller, Springer- Verlag, New York (1990); M.V. Berry and M. 

Tabor, Proc. R. Soc. London A 356, 375 (1977). 
[5] E. B. Bogomolny, B. Georgeot, M.-J. Giannoni, and C. Schmit, Phys. Rev. Lett. 69, 1477 (1992). 

[6] Aubry S. and Andr G., Ann. Israel Phys. Soc. 3 (1980) 133;M. Ya. Azbel, Zh. Eksp. Teor. Fiz. 46, 929 (1964); D. R. 
Hofstadter, Phys. Rev. B 14, 2239 (1976);D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Phys. Rev. 
Lett. 49, 405 (1982). 

[7] G. Casati, B. V. Chirikov, F. M. Izrailev and J. Ford, in Stochastic Behavior in Classical and Quantum Hamiltonian 
Systems, edited by G. Casati and J. Ford, Lecture Noted in Physics Vol. 93 (Springer, Berlin, 1979); B. V. Chirikov, F. 
M. Izrailev and D. L. Shepelyansky, Sov. Sci. Rev. 2C, 209 (1981);B. V. Chirikov, Phys. Rep. 52, 265 (1979). 

[8] S. Fishman, D. R. Grempel, and R. E. Prange, Phys. Rev. Lett. 49, 509 (1982). 

[9] I. L. Aleiner and A. I. Larkin, Phys. Rev. B 54, 14 423 (1996); Phys. Rev. E 55, R1243 (1997). 
[10] A. V. Andreev and B. L. Altshuler, Phys. Rev. Lett. 75 902 (1995). 

[11] A. V. Andreev, O. Agam, B. D. Simons, and B. L. Ahshuler, Phys. Rev. Lett. 76, 3947 (1996). 

[12] P.W. Anderson, Phys. Rev. 109, 1492 (1958). 

[13] E. Abrahams et al., Phys. Rev. Lett. 42, 673 (1979). 

[14] M.L. Mehta, 'Random Matrices', 2nd Ed., Academic Press (San Diego, 1991). 



16 



[15] J. T. Edwards and D. J. Thouless, J. Phys. C 5, 807 (1972). 

[16] V.E. Kravtsov and K.A. Muttalib, Phys. Rev. Lett. 79, 1913 (1997);S. Nishigaki, Phys. Rev. E 59, (1999) 2853. 

[17] B.I. Shklovskii, B. Shapiro, B.R. Sears,P. Lambrianides and H.B. Shore, Phys. Rev. B 47, (1993) 11487. 

[18] B.L. Altshuler, I.K. Zharekeshev, S.A. Kotochigova and B.I. Shklovskii, JETP 67 (1988) 62. 

[19] M. Schreiber and H. Grussbach, Phys. Rev. Lett. 67, 607 (1991); H. Aoki, J. Phys. C 16, L205 (1983). 

[20] F. Evers and A. D. Mirlin, Phys. Rev. Lett. 84, 3690 (2000); E. Cuevas et al., Phys. Rev. Lett. 88, 016401 (2002); A. D. 

Mirlin et al., Phys. Rev. E 54, 3221 (1996); Y. V. Fyodorov and A. D. Mirlin, Phys. Rev. Lett. 67, 2405 (1991);E. Cuevas, 

et.al., Phys. Rev. Lett. 88 016401 (2002). 
[21] M. Janssen, Int. J. Mod. Phys. B 8, 943 (1994); B. Huckestein, Rev. Mod. Phys. 67, 357 (1995); E. Cuevas, Phys. Rev. B 

68, 184206 (2003). 
[22] A. D. Mirlin, Phys. Rep. 326, 259 (2000). 

[23] F. Wegner, Z. Phys. B 36, 209 (1980); F. Wegner, Z. Phys. B 25, 327 (1976). 

[24] D. VoUhardt and P. Wolfle, Phys. Rev. Lett. 48, 699 (1982); ibtd. 45, 842 (1980). 

[25] R. Abou-Chacra, P. Anderson, and D. J. Thouless, J. Phys. C: Solid State Phys. 6, 1734 (1973). 

[26] G.M. Zaslavsky, Phys.Rep. 371, (2002) 461. 

[27] A.M. Garcia-Garcia, J. Wang, Phys. Rev. Lett. 94, 244102 (2005). 
[28] B. Hu, B. Li, J. Liu and Y. Gu, Phys. Rev. Lett. 82 (1999) 4224. 
[29] J. Liu, W. T. Cheng, C.G. Cheng, Comm. Theor. Phys. 33 15 (2000). 
[30] A.M. Garcia-Garcia, J. Wang, Phys. Rev. E 73, 036210 (2006). 
[31] A.M. Garcia-Garcia, Phys.Rev. E69 (2004) 066216. 
[32] A.M. Garcia-Garcia, J. Wang, arXiv:0707,3964 

[33] E.B.Bogomolny, C. Schmit, Phys. Rev. Lett. 93, 254102 (2004); O. Giraud, J. Marklof, and S. O'Keefe, J. Phys. A 37, 
L303 (2004). 

[34] B.L. Altshuler, L. S. Levitov, Phys. Rep. 288, 487 (1997). 
[35] D. Wintgen and H. Marxer, Phy. Rev. Lett. 60, 971 (1988). 



